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The properties of diluted Gai-^MnajAs are calculated for a wide range of Mn concentrations within 
the local spin density approximation of density functional theory. Miilliken population analyses and 
orbital-resolved densities of states show that the configuration of Mn in GaAs is compatible with 
either 3d^ or 3d^, however the occupation is not integer due to the large p-d hybridization between the 
Mn d states and the valence band of GaAs. The spin splitting of the conduction band of GaAs has a 
mean field-like linear variation with the Mn concentration and indicates ferromagnetic coupling with 
the Mn ions. In contrast the valence band is antiferromagnetically coupled with the Mn impurities 
and the spin splitting is not linearly dependent on the Mn concentration. This suggests that the 
mean field approximation breaks down in the case of Mn-doped GaAs and corrections due to multiple 
scattering must be considered. We calculate these corrections within a simple free electron model 
and find good agreement with our ab initio results if a large exchange constant {Nfi = — 4.5eV) is 
assumed. 

PACS numbers: 75.50.Pp, 75.30.Et, 71.15.Mb, 71.15.Fv 



I. INTRODUCTION 

The discovery of ferromagnetic order in diluted mag- 
netic semiconductors (DMS) made of heavily Mn-doped 
InAs and GaAs §, |, I paves the way for many new 
semiconductor spin-devices In particular the ferro- 
magnetism of Gai_a;Mna;As adds the spin degree of free- 
dom to the GaAs/(Al,Ga)As system which, in the last 
few years, has been the benchmark for new physics and 
for high speed electronic and optoelectronic devices. 

Long spin-lifetime Q and spin-coherence § in GaAs 
have already been demonstrated. Recently the feasibil- 
ity of spin-injection into GaAs using (Ga,Mn)As con- 
tacts has been proved ^ overcoming the intrinsic diffi- 
culty of injecting spins into semiconductors from mag- 
netic metals [^j. These two effects suggest that the 
GaAs/(Al,Ga)As/(Ga,Mn)As system is the best candi- 
date for injecting, storing and manipulating spins in en- 
tirely solid state devices; a valuable step towards a prac- 
tical realization of quantum computing |lO[ . 

Although there is general agreement on the carrier- 
(hole-) mediated origin of the ferromagnetism in 
(Ga,Mn)As, the detailed mechanism is still a matter of 
debate |l^, Recently Dietl et al. studied the 

ferromagnetism of III-V DMS within the Zener model, 



and obtained good agreement with existing experimental 
data using few phenomenological parameters. One of the 
key elements of the model is the mean-field Kondo-like 
coupling (j)-d Hamiltonian) between the valence band of 
the host semiconductor and the magnetic impurity 
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where N (3 is the p-d exchange constant, s is the valence 
band electron spin and S is the impurity spin. In this 
model the exchange constant, which governs the spin 
splitting of the valence band of the host semiconduc- 
tor, enters quadratically in the expression for the Curie 
temperature. Therefore its exact evaluation is crucial 
for making quantitative predictions about the ferromag- 
netism in both existing and novel materials. 

Unfortunately, in contrast with the case of II- VI 
DMS's, the experimental determination of Nfi is not 
conclusive, and both the sign and the magnitude are 
not well known, particularly for large Mn concentrations 
(x > 0.01). From the exciton splitting in the low dilution 
limit [x < 0.001) the coupling is found to be ferromag- 
netic with an exchange constant of N(3=+2.5 eV ||l4| , 
if the exchange constant for the conduction band Na is 
assumed to be -1-0.2 eV (a typical value for Mn in II- VI 
semiconductors). Refiectance magnetic circular dichro- 
ism [ p^ and magnetoabsorption experiments [T6| present 
controversial results since the absorption edge splitting 
is strongly dependent on the hole concentration (Moss- 
Burstein effect), which in turn is difficult to determine 
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from transport measurements because of a strong mag- 
netoresistance up to very high magnetic fields ||^. Mag- 
netotransport experiments are able to measure only the 
magnitude of the exchange constant and the values ob- 
tained vary from |A^/3|=3.3 eV to \NP\ = 1.5 e V M . 
Finally a recent core-level photoemission study [|19[ of 
Gao.926Mno.o74As gives Nf3=-1.2 eV if the Mn'^+ config- 
uration is assumed. Nevertheless it is worth noting that 
the raw data are compatible with both the Mn^+ and the 
Mn'^+ configurations and so one cannot make a definite 
determination of the sign of N(3. 

From a theoretical point of view, the exchange interac- 
tion for the conduction band results from a direct coulom- 
bic exchange and is expected to be ferromagnetic. In con- 
trast the exchange interaction of the valence band has a 
kinetic energy origin. It can be described as a virtual 
process in which electrons from the valence band jump 
onto the localized 3d states of Mn Therefore the 

sign and strength of the coupling depend critically on the 
population of the spin-polarized Mn d-orbitals. For less 
then half-filled d-orbitals Hund's rules suggest that the 
coupling is ferromagnetic; for half- and more than half- 
filling the coupling is antiferromagnetic. In this respect 
the case of Mn is rather peculiar. Three types of Mn 
centers in GaAs are possible. The first two can be seen 
as substitutional Mn^+ and Mn^+ respectively, with the 
former neutral (A*^ with formal Sd^ configuration) with 
spin S = 2 and the latter negatively charged {A~ with 
formal configuration 3dP) with spin S = 5/2. The third 
center is obtained when A~ weakly binds a hole, forming 
a neutral {3d^ + h) complex. The A~ center provides only 
antiferromagnetic coupling with the valence band, while 
the neutral centers can provide either ferromagnetic 
or antiferromagnetic. The ferromagnetic coupling can 
arise both from the half-filled 3d'^ center and by hopping 
through the spin polarized bound hole in the (3d^ + h) 
case pl| . 

From this brief overview it is clear that a detailed de- 
scription of the electronic structure of Mn in GaAs is 
crucial for understanding and modeling correctly the fer- 
romagnetism of Gai-^rMn^r As. In this paper we address 
this issue by calculating the ground state properties of 
(Ga,Mn)As over a range of Mn concentrations using den- 
sity functional theory (DFT) |Q in the local spin density 
approximation (LSD A). We use a numerical implemen- 
tation of DFT based oii pseudopotentials and pseudo- 
atomic orbitals ^ Although the convergence 
versus basis set with localized orbitals is more difficult 
than with plane waves (where a single parameter, the 
energy cutoff, determines the completeness of the ba- 
sis) , the method has the great advantage of being able to 
handle a large number of atoms with an accuracy com- 
parable to plane-waves methods. This allows us to in- 
vestigate various Mn dilutions without the need of large 
computer resources. Moreover the pseudo-atomic basis 



is very convenient for analysis of atomic occupation and 
orbital-resolved DOS, since no overlap integrals between 
different bases have to be calculated. 

The remainder of this paper is organized as follows. 
In the next section we provide some technical details 
about the calculation method, illustrating in particular 
how to optimize the pseudo-atomic basis set. Then we 
present our results for Gai^x^^x-^s for Mn concentra- 
tions ranging from x = 1 to x = 0.02. We analyze the 
density of states projected onto the different orbital com- 
ponents, the charge distribution around the Mn ions and 
we perform Miilliken population analyses to determine 
the occupation of the d-orbitals of Mn. In section IV we 
discuss the p-d exchange constant and we compare our 
results with that expected from the Kondo-like effective 
Hamiltonian in the mean field approximation. Then 
we illustrate how the mean field picture breaks down in 
the case of Mn in GaAs and how the LDA results can 
be explained by a simple model which includes multiple 
scattering contributions. Finally in section VI we present 
our conclusions. 



II. COMPUTATIONAL TECHNIQUE 

Since we are interested in the calculation of the elec- 
tronic properties of diluted alloys systems, we need a 
method that is able to handle with sufficient accuracy 
a large number of atoms within a periodic supercells ap- 
proach. To that purpose, we use a DFT approach based 
on pseudopotentials, and numerical localized atomic or- 
bitals as basis sets. This method, implemented in the 
code SIESTA |2^, ^ , combines accuracy and a small 
computational cost compared to other approaches with 
considerably larger computational requirements, such as 
plane waves. In this approach, however, special care must 
be devoted to the optimization of the basis set, in or- 
der to obtain the desired accuracy. In this section, we 
describe the optimization procedure used in this work. 
For all the DFT calculations presented here, we use the 
Ceperley- Alder form of the exchange-correlation po- 
tential. Self-consistency is achieved using the Pulay den- 
sity mixing scheme pTj , with a convergence criterion of 
10~^ for the change in the elements of the density matrix. 

A. Pseudopotentials 

We use the widely used scalar relativistic TrouUier- 
Martins pseudopotentials |^ with non linear core cor- 
rections |2^ and Kleinman-Bylander factorization pO| . 
The reference configurations are 4s^4p°3(i^, As^ip^Sd!^ 
and As^Ap^Sd'^ respectively for Mn, As and Ga. The 
cutoff radii for the s, p and d components of the pseu- 
dopotential are: l)Mn 2.00, 2.20 and 1.90 au, 2)As 1.90, 
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2.20 and 2.50 au and 3)Ga 2.10, 2.50 and 3.0 au. We 
check the pseudopotentials at the atomic level by com- 
paring the pseudoeigenvalues with those generated by all 
electron calculations for several atomic and ionic config- 
urations. 

In order to check the transferability of the pseudopo- 
tentials just described, we have chosen to use a plane 
waves method [Q. This allows us to perform essen- 
tially converged calculations with respect to the basis set 
(by using a sufficiently large energy cutoff for the plane 
waves). This would be difficult with the local orbitals, 
since basis set convergence is much more problematic and 
complicated in that case. In that way, we isolate pseu- 
dopotential effects from basis set effects in checking the 
pseudopotential. We have computed the equilibrium lat- 
tice constant and the band structures of both GaAs and 
MnAs, both of them showing good agreement with pre- 
viously published results js^ . 

After testing the pseudopotential, we have also cal- 
culated the band structure for a fixed localized orbitals 
basis set over a range of pseudopotentials cutoff radii. 
Our results show that, as expected, the pseudopotentials 
that yield the best results with plane waves also gives the 
best band structures with the localized atomic orbitals. 



B. Basis set: Number of (^'s 

Let us now turn our attention to the pseudoatomic ba- 
sis set. The procedure to generate the numerical atomic 
orbitals is described in Ref. The atomic orbitals are 
constructed as the product of an angular function with a 
given angular momentum I (yielding to s-type, px,py,pz- 
type, etc. orbitals), and a numerical radial function. Sev- 
eral functions with the same angular and different radial 
form can be considered to represent the same atomic 
shell, referred to as multiple-C functions. The radial 
functions are determined as follows: the first C's are 
obtained according to the scheme proposed by Sankey 
and Niklewsi Q , as the confined pseudo-atomic orbitals 
(PAO's) which result from the DFT solution of the free 
atom with the pseudopotential, and spherical potential 
of radius rc- The pseudo-wavefunction (f){r) constructed 
in such a way extends only to distances smaller than the 
cutoff radius Tc- Note that this does not correspond to a 
simple truncation, since the pseudo-wavefunction is con- 
tinuous at r = Tc. The second and successive C's are con- 
structing in the split-valence spirit. They are obtained 
by subtracting from the first C ^ function which repro- 
duces the tail of the pseudoatomic orbital for r > ruz 
and continues towards the origin as r' {a — hr"^). Here I is 
the angular momentum, a and h are parameters chosen to 
ensure the continuity, and DZ refers to "double-C" . tdz 
is chosen in such a way that the total norm beyond this 
radius has a certain value. In the present calculation we 



always fix the norm beyond tdz to 15% of the total norm, 
having noted that small variations around that value do 
not produce any significant changes in the total energy. 
Further Q's are calculated by repeating the same scheme. 
This approach is more efficient than using excited states 
of the neutral atom, which can be unbound [ p3| . 

The optimization of the PAO basis is more delicate 
than that of its plane- wave counterpart, where a single 
parameter (the energy cutoff) determines the accuracy 
and completeness of the basis. In the case of the PAO 
bases used here, several parameters determine the ac- 
curacy of the basis: the number of C's for each shell, 
the angular momentum components included, the con- 
finement radii, etc. All these must be optimized to 
achieve the required accuracy. As well as in plane- waves, 
more complete bases produce more accurate results, but 
also require larger computational resources. Since we 
are interested in describing the magnetic properties of 
(Ga,Mn)As, it is natural to choose a magnetic quantity 
as the one to monitor the convergence of our results with 
respect to the basis set quality. We study the energy 
difference ApA between the ferromagnetic and antifer- 
romagnetic alignments of two Mn atoms in a four atom 
unit cell of zincblende MnAs as a function of the basis set. 
The lattice constant is chosen to be oq — 5.8A, which is 
the critical lattice constant for the half metallic behavior 
of MnAs H. 

The first problem we address is the number of C's to 
include for each atomic orbital. We start by choosing 
double-C for the s orbitals of both Mn and As and single- 
C for the p orbitals and for the d orbitals of Mn, then pro- 
gressively increase the number of basis orbitals. The ini- 
tial cutoff radii are chosen as indicated in table || and are 
proportional to the positions of the maxima of the uncon- 
strained pseudo- wave- functions. In figure we present 



Orbital 


rc (au) 


Mn s 


6.0 


Mn p 


6.0 


Mn d 


5.0 


As s 


5.5 


As p 


5.5 



TABLE I: Initial cutoff radii used for fixing the number of ("'s. 
The radii are given in atomic unit. 

the total energies for the ferromagnetic (iSpivi) and anti- 
ferromagnetic (-Eaf) alignments, and ApA, for the differ- 
ent PAO bases listed in table |ll[ From the picture two 
important conclusions can be reached. First, according 
to the usual variational principle the total energies for 
both ferromagnetic and antiferromagnetic configurations 
decrease with enlarging the basis. Second the split be- 
tween the ferromagnetic and antiferromagnetic configu- 
rations is significantly reduced by using triple-^'s for the 



4 



Basis 


Nc(As s) 


Nc(As p) 


Nj(Mn s) 


Nc(Mn p) 


Nj(Mn d) 


1 


2 


1 


2 


1 


1 


2 


2 


1 


2 


1 


2 


3 


2 


1 


2 


2 


2 


4 


2 


2 


2 


2 


2 


5 


2 


2 


2 


2 


3 


6 


2 


2 


3 


2 


3 


7 


2 


3 


2 


2 


3 


8 


2 


2 


2 


3 


3 



TABLE II: Summary of the bases used in figure ^ In the first 
column the indicator of the basis, the following columns show 
the number of ^ (N^) for each orbital. 



d orbitals of Mn (note the large decrease of ApA when 
going from basis 1 to basis 2 and from basis 4 to basis 5) 
and doublc-C for the p orbitals of As. This sensitivity of 
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FIG. 1: (a) Efm, Eaf and (b) Afa for the basis of table |l[ 
Note the decreasing of the total energies as the basis becomes 
more complete. Afa saturates for the basis 5. 



the magnetic phase stability to the As-p and Mn-d basis is 
consistent with the magnetism in MnAs being driven by 
strong p-d hybridization. Since we are mainly interested 
in the magnetic properties of diluted systems describable 
by very large unit cells we decide to use double-^ for all 
the orbitals except the d orbitals of Mn for which we use 
triple-C (this is the fifth basis set in table 2). Note that 
we can afford to use triple-^ for Mn-d since few Mn ions 
are present in the cell. In contrast the use of larger basis 
sets for Ga and As yields a more dramatic increase of the 
size of the computations. 



C. Basis set: Cutoff radii 

Next we turn our attention to the choice of the cut- 
off radii of the basis sets. For free atoms, the optimum 
cutoff radius of any orbital (as the one which minimizes 
the energy) is infinite, since that case corresponds to no 
confinement potential, which yields to exponential tails 
for all the atomic wavefunctions. However, in solids this 
criterion does not hold, since the lack of vacuum and 
the presence of the crystal potential tend to confine the 
atomic wavefunctions more than in the free atom. In this 
situation, the confinement of each PAO should be opti- 
mized to minimize the total energy. This procedure has 
shown in other systems like bulk bcc Fe that a fi- 
nite and relatively small confinement radius can provide 
lower energies and therefore more accurate bases than 
long values of rc- These calculations also show that the 
optimum confinement radius of each PAO depend very 
much on the particular orbital. In our case, however, it 
would be too complex to optimize the rc for all the or- 
bitals in our system, due to the large number of them. 
Instead, we have followed a simpler approach. We vary 
the cutoff radii of table | uniformly by multiplying all the 
radii by a common scaling factor t. A somewhat similar 
criterion is to use as variational parameter the orbital en- 
ergy shift Aii^pAo, which is the energy increase that each 
orbital experiences when confined to a finite sphere and 
can be used as single parameter to test the convergence 

In figure H we present E'fm, ^'AF and Afa as a function 
of t for the fifth basis set of table 2. The last point 




FIG. 2: (a) ii^FM, Eaf and (b) Afa for the fifth basis set of 
table ^ as a function of the scaling parameter t. Note that 
Afa saturates sooner than the total energies for ferromagnetic 
and antiferromagnetic alignment. 

in figure || (t = 1.41) corresponds to a basis with an 
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orbital energy shift Ai?pAo of O.OOlRy. A convergence 
of O.OOlRy has been successfully used to describe the 
magnetic properties of Ni clusters on Ag surfaces 
and is considered an optimal value for the convergence. 
However in our case we prefer to use smaller cutoff radii 
in order to reduce the computation time. From figure |^ 
it is clear that the saturation of Ap^ occurs for shorter 
radii than those required to converge the total energies. 
We therefore decide to fix the cutoff radii to i = 1 noting 
that E-pM and Eaf differ from the value obtained for 
Ai;pAo=0.001Ry by only 0.04% and that Afa differs by 
only 2%. 



D. Comparison with previous calculations 

We further test our basis set by computing the energy 
split between the Mn d states with e and t2 symmetry at 
the r point, and the dependence of the magnetization on 
the lattice spacing for zincblende MnAs. These two tests 
give an indication of the accuracy of the p-d exchange, 
which is a dominant interaction in (Ga,Mn)As. In fact 
at the r point the t2 states are coupled with the As-p 
states while the e are decoupled, and their splitting is 
determined by the p-d coupling. 

In figure ^ we present the energy split Ai?e-t2 = 
Ee ~ Et2 as a function of t for both the spin directions. 
The 6-^2 split converges monotonically and there is a 




FIG. 3: e-t2 energy split at the F point for MnAs; (a) majority 
spin, (b) minority spin. 

variation of only ^2% going from t — 1 to t — lAl 
(A£'pAo=0.001Ry). If we now compare this result with 
our previously published results obtained with plane- 
waves we notice that our present results give an e-t2 split- 
ting around 50meV less than the plane-wave splitting for 
both spins. This is roughly the same discrepancy found 



for ApA- As we have just shown, such a deviation from 
the plane-wave calculation cannot be lifted by increasing 
the size of the basis, since this does not produce vari- 
ations larger than 2%. A possible origin of such a dis- 
agreement may be the slightly different pseudopotentials 
used. 

In figure ^ we present the magnetization as a function 
of the lattice constant for MnAs and we compare it with 
that obtained previously in our plane-wave calculations 
P4. The agreement is quite good; the transition to the 
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FIG. 4: Magnetization as a function of the lattice constant 
for zincblende MnAs. 

half-metallic state is correctly predicted for ao=5.8A and 
also the dependence of the magnetization on the lattice 
spacing is well reproduced for ao > 5.3A. For smaller 
lattice spacings the two calculations disagree with a ten- 
dency of the PAO basis to over-stabilize the ferromag- 
netic phase. This is a quite general behavior that we 
found also with other basis sets. In particular for bases 
with a small number of C's this effect is magnified. For 
instance the magnetization calculated with the first basis 
set of table || for ao = 5.0A is 2.8/j,b, while that calcu- 
lated with our optimized basis is ^ 1.4/j,b. However we do 
not believe that our plane- wave and localized basis calcu- 
lations should necessarily agree for strongly compressed 
unit cells since the portability of the two pseudopoten- 
tials used will be different. 

Finally we check the ability of our optimized basis set 
to describe the electronic and structural properties of 
GaAs, which forms the matrix where the Mn ions are 
included in (Ga,Mn)As. We find an equilibrium lattice 
constant of ao=5.635A which is remarkably close to the 
experimental one. Moreover the bandstructure is very ac- 
curate; a comparison of our calculated eigenvalues with 



existing calculations is presented in table [II. 

In summary, we are confident that the results which 
we obtain using the numerical atomic orbitals method 
with the combination of TrouUier-Martins pseudopoten- 
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tials and the basis set described in this section, are in 
good agreement with LSDA results obtained using other 
techniques. 



III. ELECTRONIC CONFIGURATION OF 
GaAsrMn 

In this section we study the electronic structure of Mn 
in GaAs. We consider large GaAs cells (from 2 to 96 
atoms) in which we replace one Ga atom with a Mn atom 
(Mnca)- We use 18 k points in the corresponding irre- 
ducible Brillouin zones for all the supercells and over 1000 
k points for the primitive zincblende unit cell (2 atoms). 
Since the cell contains only one Mn atom and we use peri- 
odic boundary conditions the Mn atoms are forced to be 
ferromagnetically aligned. For the smaller cells (32 and 
48 atoms) we perform several simulations changing the 
shape of the unit cell. This is equivalent to investigating 
different arrangements of the Mn atoms with respect to 
each other. We find that, although the general proper- 
ties do not change, different Mn ion arrangements in the 
cell result in different total energies. For all the calcula- 
tions we assume the GaAs experimental lattice spacing 
ao=5.65A. 



A. Partial DOS and charge density distribution 

We start by analyzing the orbital-resolved density of 
states. In figure |5| we present as an example the DOS 
obtained for a 64 atom unit cell with one Mnga sub- 
stitution. Similar features are obtained for both higher 
and lower Mn concentrations. Far from the Fermi en- 
ergy the DOS remains close to the DOS of GaAs (see 
figure ^), with a lower energy As-s band and the Ga(s)- 
As{p) valence band. At the Fermi energy the situation is 
markedly different. The majority spin band has a rather 
sharp peak, characteristic of a narrow band, while the 
minority spin has a gap. Such a band structure is the 
signature of a half-metallic material. The total magne- 
tization of the cell is 4//b- Our calculations for higher 
and lower Mn concentrations show that the magnetiza- 
tion does not change with the Mn concentration. In the 
minority band the corresponding peak is shifted to higher 
energy and is very close to the edge of the GaAs conduc- 
tion band. If we now consider the DOS projected onto 
the different orbital components of the wave- function and 
we look at the e and t2 d-states of Mn some interesting 
features appear. The majority band exhibits two broad 
peaks between -4eV and -leV below the Fermi energy 
with strong e and t2 component respectively. In addition 
there is a rather narrow t2 peak at the Fermi energy. In 
contrast the minority band has almost no c?-character be- 
low but instead has two sharp e and t2 peaks around 
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FIG. 5: Partial density of state for Gai_a;Mni,As for x = 0.3 
(one Mnca in a 64 atom GaAs cell): (a) majority spin, (b) mi- 
nority spin. The vertical line denotes the position of the Fermi 
energy. The dashed and dotted lines represent respectively 
the projection of the DOS onto the Mn-d t2 and e orbitals. 




E(eV) 

FIG. 6: DOS for GaAs. 

leV above Ep. The different peak widths reflect the dif- 
ferent degrees of hybridization of the Mn-d band with 
the GaAs bands. The hybridization is much stronger for 
states far below the Fermi energy. 

In order to have a better understanding of the p-d hy- 
bridization in diluted (Ga,Mn)As, we present in figure 
1^ the evolution of the Mn-d peaks as a function of the 
Mn concentration. The most relevant feature is that for 
the sharp peaks in both the majority and minority bands 
(columns (b) and (c)) the relative intensity of the d com- 
ponent of the DOS is independent of the Mn concen- 
tration. Therefore those portions of the DOS must be 
derived almost entirely from the Mn impurity and its 
four neighboring As atoms. This can also be seen by 
looking at the DOS projected onto the p states of the 
four atoms tetrahedrally coordinated with Mn (the dot- 
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FIG. 7: Total and orbital-resolved DOS for different Mn con- 
centrations. The three columns correspond to different en- 
ergy regions and spins: (a) majority band between -4 and 
-leV below E-p (broad Mn-d peaks), (b) majority band at 
the Fermi energy (sharp Mn-d peak), (c) minority band leV 
above -Ep. The four rows indicate different Mn concentra- 
tions: (1) a;=0.06 (1 Mn in 32 atoms), (2) 2;=0.04 (1 Mn in 
48 atoms), (3) a;=0.03 (1 Mn in 64 atom) and (4) 1=0.02 (1 
Mn in 96 atoms) . The solid lines denote the total DOS while 
the dashed, dotted and dot-dashed lines denote the DOS com- 
ing respectively from the t2 states of Mn, the e states of Mn, 
and the As(p) states of the four As atoms neighboring the 
Mn impurity. Note that the states of columns (b) and (c) 
do not scale with concentration, indicating strong local hy- 
bridization. 



dashed line of figure 0) . In summary our analysis shows 
that the MnAs4 complex accounts for most of the DOS 
at the valence band edge for the majority band and at 
the conduction band edge for the minority. 

In contrast the Mn-d states far below E'p results from 
strong coupling with the p orbitals of all the As atoms 
of the GaAs cell. This can be clearly seen in figure |[ 
where we present the charge density isosurface plots cor- 
responding to the three DOS of figure |^. Figure || shows 
that the charge corresponding to states at the edge of the 
GaAs band gap is localized around the MnAs4 complex 
(figures ||b and ||;), while the remaining Mn-d states are 
hybridized with all the As-p orbitals (figure |^a) . 

We now turn our attention to the distribution of the 
magnetization around the Mn ion. The magnetization 
around one atom placed at Ro is calculated as 



M(R,)= f [pT(r-Ro)-px(r-Ro)]dr, 



(2) 



where flu. is a sphere of radius Ri and is the charge 
density for the spin a. The charge density is calculated on 
a real-space grid by evaluating the localized orbitals on 












c 













FIG. 8: Charge density isosurfaces for three states shown in 
figure ^. In this case we consider only a;=0.06 (1 Mn in a 
cubic 32 atom cell), (a) Majority band between -4 and -leV 
below Ey (broad Mn-d peaks), (b) majority band at the Fermi 
energy (sharp Mn-d peak), (c) minority band leV above Ey- 
The Mn ion is in the center of the cell. 



such a grid gj. Of course M{Ri) depends on the cutoff 
radius Ri. In figure |9| we present the magnetization of 
Mn, and of the first and second As nearest neighbors of 
Mn, as a function of Ri. 
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FIG. 9: Magnetization profile as a function of the integration 
radius Ri for (Ga,Mn)As with a;=0.3 (1 Mn ion in a cubic 64 
atom GaAs cell), (a) Mn, (b) first nearest As atom to Mn, 
(c) second nearest As atom to Mn. 
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In the case of Mn the magnetization saturates 
for i?i=4.0 au and remains almost constant up to 
Ri ~18.0 au when the next Mn shell is encountered. 
Hence we can easily deduce that the Mn magnetization 
is 4 /zb which is the saturation value. In contrast the 
magnetization around the As ions shows a negative min- 
imum (between 2 and 3 au from the As ion, depending 
on the position of the As ion relative to the Mn ion) fol- 
lowed by a sharp increase. The minimum corresponds 
to a negative spin polarization with respect to the Mn 
and the following magnetization increase occurs at dis- 
tances where the polarization of the neighboring atoms 
starts to be included in the integration. In the case of 
the first nearest neighbor this magnetization increase is 
due mainly to the spin polarization of Mn (figure ^) 
and in the case of second nearest neighbors it is due to 
the four Ga ions coordinated to As (figure ||c). It is in- 
teresting to note that the polarization of the As ion is 
always negative with respect to that of Mn. This means 
that Mn and As are antiferromagnetically coupled. The 
values of the spin polarizations of As at the minima are 
—0.03 /iB and —0.005 fis respectively for first and sec- 
ond nearest neighbors. These values of polarization are 
similar to those already published for GaAs/MnAs su- 
perlattices calculated with a first principle LMTO-ASA 
method [41]. It is worth noting that we did not find 
any sizeable changes in the magnetization per atom as a 
function of the Mn concentration for all the concentra- 
tions studied. 




E(eV) 

FIG. 10: Orbital resolved DOS for zincblende MnAs with lat- 
tice spacing ao = 5.65A. The vertical line denotes the position 
of the Fermi energy. 

Finally we compare the orbital resolved DOS of 
(Ga,Mn)As with that of zincblende structure MnAs. In 
figure we present the DOS for zincblende MnAs with 
the lattice spacing of GaAs (5.65A), which is the same 
lattice spacing that we used for diluted (Ga,Mn)As. For 



this lattice spacing MnAs is not an half-metal since the 
Fermi energy in the minority band cuts through the con- 
duction band edge, mainly dominated by d electrons. Al- 
though the total DOS is different the projection onto the 
d-orbitals closely resembles that of diluted (Ga,Mn)As 
(see figure In particular there is a large occupation of 
the d-orbitals in the majority band, while in the minor- 
ity band only the bands of t2 symmetry are occupied as 
a result of the hybridization with the As-p states at the 
edge of the valence band. The magnetization integrated 
around Mn ions is smaller in zincblende MnAs at this 
lattice constant than in (Ga,Mn)As. For a lattice spac- 
ing of ao=5.65 A we find a Mn polarization of 3.79 /ib, 
compared with 4.0 /zb of (Ga,Mn)As. In contrast the po- 
larization of As in zincblende MnAs is considerably larger 
with an integrated magnetization of about —0.17 /^b- We 
also note that on increasing the lattice spacing the spin 
polarization of Mn increases, but the polarization of the 
As is largely unchanged. For instance at ao=5.80 A we 
find 4.04 /ib and —0.17 /ib respectively for the Mn and 
As magnetizations. This suggests that the polarization 
of Mn is related to the ionicity of the bond with As. 

A more quantitative comparison of the zincblende 
MnAs with diluted (Ga,Mn)As can be obtained by per- 
forming Miilliken population analyses We de- 
scribe the results of such analyses in the next section. 



B. Miilliken population analysis 

We perform Miilliken population analyses Q in 
order to compare quantitatively the orbital occupations 
of (Ga,Mn)As at different dilutions. The Miilliken pop- 
ulation analysis is a convenient way to separate dif- 
ferent contributions to the total charge density. Sup- 
pose we have a system described by the wave-function 
(j> — ciipi + C2?/'2, where ipa is a localized function and 
Cq is the corresponding amplitude. Then if the state (f> 
is occupied by N electrons, the total occupation can be 
written as 



N ^Ncl + 2NC1C2S12 + Nc 



2 ' 



(3) 



where 5*12 is the overlap integral, J ipi'ip2dv. Miilliken 
defined the sub-populations Ncf and Nc2 as net popu- 
lations and 2NC1C2S12 as overlap population. Moreover 
if the overlap population is equally splitted between the 
two wavefunctions we obtain respectively Nc1+NciC2Si2 
and Nc2 + NC1C2S12, which are referred as gross popula- 
tions. In what follows we always refer to the gross popula- 
tion. If the functions ipa represent orbital components of 
the angular momentum, then the populations correspond 
to orbital populations and the overlap population is the 
orbital overlap population. Similarly if the functions 
ipa are the atomic wave functions for the atom a then 
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we obtain the atomic populations and the atomic over- 
lap population. We also define Miilliken atomic charge 
as the diiTerence between the gross atomic charge (ie 
eNc2 + eNciC2Si2 with e the electronic charge) and the 
valence charge of the isolated atom. It is widely accepted 
that the absolute magnitude of the atomic charges can 
depend strongly on the basis set in which they are cal- 
culated ^J] . However relative values of Miilliken popu- 
lations can provide useful information when comparing 
different systems (for instance the amount of covalency 
in semiconductors) ]4^ . 

We start the analysis by calculating the Miilliken 
atomic charges for Ga, Mn and As in GaAs, MnAs and 
(Ga,Mn)As at different concentrations (table IV). In the 



case of (Ga,Mn)As for each atomic species we present the 
average values over the cell. Since we have already shown 
that the four As atoms coordinated with Mn (which we 
denote by As-^^) have quite different properties than the 
remaining As atoms, we calculate their average atomic 
charge separately. The table shows clearly the local char- 
acter of the MnAs4 center. We note that the average 
Miilliken charges of Ga and As closely resemble those of 
GaAs, particularly for low Mn concentrations. Of course 
in the extremely diluted limit GaAs:Mn one expects the 
average charges of Ga and As to be exactly those of 
GaAs. In contrast, the average Miilliken charge of the 
four As atoms coordinated with the Mn impurity does 
not change with concentration, confirming that the elec- 
tronic structure of the MnAs4 complex is not affected by 
the concentration. It is also interesting to note that these 
As atoms have small positive atomic charges whereas the 
other As atoms have negative atomic charge. A positive 
As atomic charge is also found in zincblende MnAs, al- 
though in that case its magnitude is much larger. The 
transition from GaAsiMn to zincblende MnAs with in- 
creasing Mn concentration is reflected in the increase of 
Miilliken charge on the Mn atoms. Therefore the Mn- 
As bond becomes more ionic when the Mn concentra- 
tion is increased. This picture, together with the almost 
complete occupation of the d shells in the majority band 
discussed in the previous section, is consistent with mod- 
eling Mn in GaAs as an impurity center composed of 
a negatively charged Mn ion in a conflguration, and a 
weakly bound hole {d'' + h) |46[ |47|. The increase of Mn 
concentration, and the consequent increase of hole con- 
centration, reduce the binding energy of the bound hole 
due to the partial screening of the potential. Therefore an 
increase of the Miilliken charge of Mn with concentration 
is expected. Nevertheless the agreement is only qualita- 
tive and a deflnitive prediction based solely on Miilliken 
analysis cannot be made. 

We now turn our attention to the orbital population. 
In table ^ we present the orbital populations for the p 
orbitals of As and the d orbitals of Mn in MnAs and 
(Ga,Mn)As for both spin orientations. As before, we 



distinguish between the As''^^ atoms and the remaining 
As atoms. We do not report the orbital populations for 
Ga, for the s orbital of As and for the s and p orbitals of 
Mn since they are not relevant for the discussion. 

Several important aspects can be pointed out from the 
table. The total population for the d orbital of Mn is 
around 5.5 electronic charges for all the systems studied. 
We do not expect integer values for the orbital popula- 
tion since strong p-d hybridization is present. The to- 
tal overlap population for zincblende MnAs is about 0.6 
electronic charges and this can be considered to be the 
uncertainty on the determination of the orbital popula- 
tion. This give a Mn d orbital occupation of 4.6 ± 0.6 
and 0.7 ± 0.6 for the majority and minority bands re- 
spectively. Although the orbital population is not an ob- 
servable quantity and its absolute value may be affected 
by the choice of the basis set, we can conclude that the 
atomic configuration of Mn in GaAs is compatible with 
both and 3d^ . This is in agreement with recent x-ray 
absorption magnetic circular dichroism experiments p8| , 
where the data are interpreted by assuming a Mn config- 
uration consisting of 80% Mn 3d^ and 20% Mn M^. It 
is interesting to note that by decreasing the Mn concen- 
tration there is an increase of the polarization of the d 
orbital of Mn (the orbital population is enhanced in the 
majority band and reduced in the minority). This seems 
to be in favor of the A~ configuration in the limit of 
high dilution, as is reported extensively in the literature 

MM- 

also shows clearly that there is antiferromag- 
netic coupling between the Mn d and As p orbitals. The 
As p orbitals in fact have quite a large spin polarization 
opposite to that of Mn. This cannot be due to the overlap 
components of the orbital population, which would give 
the same polarization as that of Mn. It is also interesting 
to note that the spin polarization is much larger among 
the As^^ atoms for which it is almost insensitive to Mn 
concentration, than among the other As atoms. As ex- 
pected it is still smaller than the As spin polarization 
in zincblende MnAs. Nevertheless we also see that the 
other As atoms have a small antiferromagnetic polariza- 
tion of the p orbital, which decreases with concentration 
as expected. This is in very good agreement with the 
magnetization data presented in the previous section. 



IV. EXCHANGE COUPLING 

As we pointed out in the introduction, the evaluation 
of the exchange constant N(3 is crucial for predicting the 
thermodynamic properties of (Ga,Mn)As. In this sec- 
tion we provide a theoretical estimate of the exchange 
constant and study its dependence on the Mn concentra- 
tion. We begin by briefly describing the effect of the sp-d 
exchange on the bandstructure of the host semiconductor 
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in the mean field approximation. Our starting point is 
the commonly used sp-d exchange Hamiltonian pO| 



H 



sp- 



n,k,k' 



E 



(4) 



where J^p ''(k, k') is the exchange integral of the band 
electrons (n, k) and (n, k') with the Mn local spin Si, c|jj^ 
and c„k are the creation and annihilation operators for an 
electron in band n with Bloch vector k. The sum extends 
over the valence (n = v) and conduction [n = c) bands 
of GaAs, and all the localized spins labelled by the index 
i. If we neglect interband terms which are negligible and 
replace the spin Si by the average spin (S) proportional 
to the magnetization we restore the translational invari- 
ance of the system. Therefore the expression (||) becomes 
diagonal in k and can be written as a function of the Mn 
fraction x and the cation concentration N as 



H. 



1 



sp—d 



(5) 

with t (j) indicating the up spin (down spin) direction 
with respect to the mean field spin {S) and J^^'^{k) = 
J^P~'^{'k, k). If we now restrict our analysis to the band- 
edge (r point) and define a = J^p-'^{0) and /3 = J^P-'^iO) 
we obtain the equations 



sp—d 



'^xN{S)a{clQ^Cco^ 



,t 



(6) 



We calculate the exchange constants directly from 
the conduction band-edge (valence band-edge) spin- 
sphttings AE'= ^ El - E^^ (Ai;" ^ El - E^) as follows 

Na^ AEyx{S) , N(3 = AEyx{S), (8) 

where (S) is half of the computed magnetization per Mn 
ion. In order to evaluate the parameters in equations (||) 
we compute the band structure around the T point for 
large GaAs cells with a single Mn impurity. In figure 
[lT| we present as an example the results for a cubic cell 
containing 64 atoms. Since we are mainly interested in 
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FIG. 11: Band structure and orbital resolved DOS at the F 
point for Gai_2,Mna;As with x=Q.3 (1 Mn ion in a cubic 64 
atom GaAs cell): (a) majority band, (b) minority. 



for the conduction and valence bands respectively. We 
note that the same analysis can be carried out by assum- 
ing that the Mn impurities form a perfect ferromagnetic 
crystal. In such a case the derivation of equations (||) is 
identical to that given here if the magnetic moment per 
Mn atom is used for the mean field spin (S). Equation (||) 
relates the spin-splitting of the conduction and valence 
bands to the exchange integral calculated at fc = 0. This 
quantity is usually extracted in optical magnctoabsorp- 
tion experiments from the spin-splitting of the cxciton 
lines. For instance the Zeeman splitting of the heavy 
hole exciton transition Ei is 

Ei=x{S)N{(3-a). (7) 

Other transitions give different combinations of a and 
P, which can then be determined. Note finally that the 
spin-splitting of both the valence and conduction bands 
in the mean field approximation is linear with the Mn 
concentration x. 



the r point we consider the band structure only along the 
direction (^-,0,0) ^ (0,0,0) ih-A-A-) with 

^ 8 CO ' ' ' ^ ' ' ' ^ 8 CO ' 8 CO ' 8 CO ' 

Co the unit vector of the cubic cell. We indicate these 
two directions respectively as X and M. In figure |l^ we 
also plot the orbital resolved DOS at the F point. This 
shows clearly that the valence band edge has mainly As-p 
character with contributions also from the t2 Mn-d states 
due to hybridization, while the conduction band edge is 
formed by Ga-s states. In this way the spin-splitting is 
easily computed. 

We consider different Mn concentrations and for the 
smaller unit cells (larger concentrations), different geo- 
metrical arrangements. We find that the spin-splittings 
of both the conduction and the valence bands are depen- 
dent on the relative positions between the Mn ions, with 
variations of up to 20%. In particular for the same Mn 
concentration we find large splittings when the Mn ions 
are clustered, and smaller splittings for homogeneously 
diluted systems. More details on the dependence of the 
exchange constant on the spatial arrangement of the Mn 
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ions will be published elsewhere ||5^. In the following 
we consider only cells which maximize the separation be- 
tween the Mn ions (uniform Mn distribution). 

In figure |lj we present the spin-splitting for the con- 
duction and for the valence band as a function of the Mn 



concentration and in table VI we list the corresponding 
exchange constants. Consider first the conduction band. 



> 




> 




FIG. 12: Spin splitting of the conduction (a) and valence 
(b) bands as a function of the Mn concentration x for 
Gai-xMna;As. The corresponding exchange constants are 
listed in table 0. 

Although the spin splitting shows large fluctuations with 
X, there is no systematic variation with the Mn concen- 
tration. With the caveat that DFT is a ground state the- 
ory and therefore does not describe accurately the con- 
duction band, from table VI one can conclude that the 



coupling {s-d coupling) between the conduction band of 
GaAs and the Mn impurity is ferromagnetic. Also it is 
independent of predicted by mean field theory, and 
has the value Na ^0.2 eV. Note that ferromagnetic cou- 
pling is expected since in the case of the conduction band 
the only exchange is direct and also that the value of the 
exchange constant Na is very close to that usually found 
in II- VI semiconductors [ p3[ . 

The situation is quite different for the valence band. 
First of all we see that the spin-splitting of the valence 
band edge is much larger than the typical absorption edge 
splitting in magnetoptical experiments [ p^ . For instance 
if we compare the results for x = 0.032 of reference [^g) 
with those of figure 
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for the same concentration, we 
find that our calculated value is about four times larger 
than that obtained experimentally. However it is impor- 
tant to point out that in our calculation all the Mn ions 
contribute to the ferromagnetism. In contrast in real 
systems only a fraction of the Mn ions are ferromagneti- 
cally aligned, and the typical magnetization curves have 
a large paramagnetic component which does not saturate 
even at very high magnetic fields . This has been con- 
firmed by recent x-ray magnetic dichroism measurements 
Assuming a mean field picture, this suggests that 
the mean field spin calculated here is much larger than 
that present in actual samples. Turning the argument 
around we can conclude that our results are consistent 



with experiments if we assume that in the latter the ef- 
fective Mn concentration (contributing to the ferromag- 
netism) is only 1 /4 of the nominal concentration. 

A second important point is that, although figure |l2| b 
seems to suggest a linear behavior of the spin-splitting 
AE'^ with X according to equation (||) , a closer look at ta- 
ble VI reveals that the exchange constant Nf] is strongly 
dependent on x. Specifically, Nf3 increases with decreas- 
ing Mn concentration, a behavior already well known to 
occur in Cdi^ajMn^^S ||5^, |6[ This dependence of 

the exchange constant on x could be due to two possible 
reasons: i)the actual p-d coupling is not Kondo-like, or 
ii)the mean-field approximation that leads to equation 
(P) is not valid. Blinowski and Kacman studied the ki- 
netic exchange interaction of various 3d metal impurities 
in zincblende semiconductors p8| . By applying canoni- 
cal transformations to the p-d hybridization Hamiltonian 
|59|| they evaluated the effective exchange interaction be- 
tween the valence band and the d states of the impu- 
rity. They found that for the and configurations 
of the impurity the effective exchange has a Kondo-like 
form, while there are other non Kondo-like contributions 
for the 3c?'* case. From the Miilliken analysis we can 
rule out this latter configuration, and conclude that the 
effective exchange is indeed Kondo-like. Therefore the 
dependence of iV/3 on x is suggestive of the breakdown 
of the mean-field approximation. 



V. BREAKDOWN OF MEAN FIELD 
APPROXIMATION 

The main hypothesis sustaining the mean field approx- 
imation is that the potential introduced by the Mn ions 
is weak with respect to the relevant band-width. This 
seems to be true in most of the II- VI semiconductors, 
however in the case of (Cd,Mn)S such a hypothesis breaks 
down and an apparent strong dependence of the exchange 
constant on the Mn concentration is found |5^, pTf . 
The case of Mn in GaAs looks very similar. We recall 
that for the very diluted limit there is some evidence of 
the Mn ion being able to bind a polarized hole |^ . This 
suggests that the potential created by Mn in GaAs may 
be strong and hence the mean field approximation breaks 
down. 

Benoit a la Guillaume et al. |Q calculated the correc- 
tions to the mean field approximation using a free elec- 
tron model, with the magnetic impurities described by 
square potentials. The energy was calculated within the 
Wigner-Seitz approach which is applicable only to the 
case of perfectly periodic crystal. Although the model 
has been refined ^ the main findings are still valid. 
Here we illustrate briefly the model and we use it for 
computing the exchange constant. 

We consider a free electron model with effective mass 
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m*, and uniformly distributed magnetic impurities de- 
scribed by the potential 



U{r) = W{r) - J{r)S ■ s. 



(9) 



Here W{r) is the spin independent substitutional po- 
tential and J{r) is the p-d coupling between the free 
electron spin s and the impurity spin S. We further 
assume that J(r) and W{r) have the same square po- 
tential shape, and that all the impurity spins are ferro- 
magnetically aligned. This leads to U{r) = Uq 6{r — h) 
and also = {W T 5/AJ)0{r - b), when S = 

5/2 is considered. Finally the energy is calculated by 
solving a transcendental equation obtained by impos- 
ing the appropriate boundary conditions We de- 
fine 6{x,Uo) = E{x,Uo)/Eynfa{x,Uo) as the deviation of 
the computed energy E{x, Uq) from the mean field en- 
ergy E^faix, Uo) = VNx where V = ^b^Uo and Nx 
is the Mn density. In figure |l^ we present 5{x) as a 
function of x for different potentials 77 = Uo/\Uc\, where 
Uc = —{TTh/2by /2m* is the binding potential. We no- 
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FIG. 13: Dependence of the correction factor 5{x) on the 
impurity concentration x for a range of potentials. Note that 
the deviation from mean field theory is larger for strongly 
attractive potentials. 

tice that the corrections to mean field theory are large 
for small x and decrease with increasing x. In particu- 
lar the mean field approximation breaks down when the 
potential is attractive and close to the binding potential 
(rj — > 1), while it is reasonably good for repulsive poten- 
tials. We also note that the mean field approximation 
is recovered in both the limit of large Mn concentrations 
{x — > 1) and in the limit of weak potential (77 — > 0). 
This general behavior can qualitatively explain our LDA 
results. Consider in fact the band edge spin-splitting 

AE^ix) = [{W + 5/4J) S{x, r]^) - 



{W -5/4.7) Six, r]'')]Nx , (10) 

where rj^'^^^ = (IF T 5/4J)/|C/c|. By comparing equation 
( p^ ) with figure |l^ one can see that for small x the spin- 
splitting is largely enhanced with respect to its mean field 
value. The deviation diminishes on increasing x, and van- 
ishes in the limit of complete Mn substitution [x — 1). 
Note that the application of the mean field approxima- 
tion at every x gives rise to an apparent increase of the ex- 
change constant with the Mn concentration. This agrees 
with our LDA results. 

It is also worth noting that the deviation from mean 
field theory is larger if the spin asymmetry of the poten- 
tial Uq is large. In particular the spin-splitting is largest 
when the potential is attractive for one spin species and 
repulsive for the other. In the opposite limit, when the 
mean field approximation is valid (6 — > 1), the equation 
( [l0| ) reduces to the usual expression 



AE^ix) = -N(3x , 



(11) 



where we have defined the exchange constant N (3 = 
iV(47r&3/3)J. 

In order to compare with experiments we perform a 
fit of our LDA data. We consider &, W and J as fitting 
parameters, with b varying between the cation-anion and 
the cation-cation distance and W and J chosen so that 
no bound holes are present. This last restriction takes 
into account the lack of any experimental evidence for 
bound holes at the concentrations investigated here. The 
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FIG. 14: Band edge spin-splitting: the circles represent our 
LDA data and the straight line the fit obtained with the 
model discussed in the text for the parameters h = 3. 6 A 
J = -1.05eV {NP = -4.5eV) and W = -0.027eV. 

fitting procedure yields b — 3.6A and N(3 — — 4.5eV, 
although values in the range 3.6A < b < 3.9A and 
-4.9eV < N(3 < -4.4eV fit equally well. It is important 
to note that for all the parameter sets which give a good 
fit one spin hole is nearly bound while the other feels a 
weak repulsive potential. Our best fit is presented in fig- 
ure 14 , Despite the roughness of the model the agreement 
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is reasonably good. It is interesting to point out that the 
model seems to underestimate the spin-splitting for small 
X and overestimate it for large x. This is not surprising; 
in the model we assume that the potential induced by the 
magnetic impurity does not depend on the impurity con- 
centration. This is in general true for Mn in II- VI semi- 
conductors, where Mn provides only a local spin. In the 
case of the III-V semiconductors, however, Mn acts both 
as an acceptor and as source of localized spins. There- 
fore it is natural to expect a progressive screening of the 
Mn potential with concentration due to the increase of 
the hole density. This effect, which is responsible for the 
lack of bound holes in low diluted (Ga,Mn)As, further 
reduces the deviation from the mean field approximation 
for large x, and better agreement with our LDA data 
may be found. 

VI. CONCLUSIONS 

We have investigated theoretically the magnetic prop- 
erties of Gai-ajMuj^As with dilutions ranging from x=l 
to X—Q.Q2. We found that Mn in GaAs has an atomic 
configuration compatible with both and 3d^, and that 
the total occupation is not integer because of the strong 
p-d coupling with the valence band of GaAs. Such a 
coupling is antiferromagnetic with a remarkably large 



exchange constant. We have shown that the exchange 
constant has an apparent dependence on the Mn concen- 
tration. This suggests that the generally used mean field 
approximation breaks down, since the potential induced 
by the Mn ions in GaAs cannot be treated perturba- 
tively. Using a simple free-electron model we have calcu- 
lated the corrections to the mean field expression for the 
spin-splitting of the GaAs valence band and found a good 
agreement with the LDA calculations. Further study is 
needed to determine the dependence of the spin-splitting 
on the confinement of the Mn ions in the case of highly 
ordered alloys. 
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GaAs 


Fi 


Fi 


Xi 


Xs 


X3 


Xi 


L2 


Li 


La 


Li 


PAOcxp 


-12.91 


0.57 


-10.53 


-7.05 


-2.84 


1.71 


-11.25 


-6.86 


-1.26 


1.34 


PAOthc 


-12.99 


0.66 


-10.56 


-7.10 


-2.88 


1.88 


-11.30 


-6.92 


-1.27 


1.39 


LAPW 


-12.80 


0.29 


-10.29 


-6.89 


-2.69 


1.35 


-11.03 


-6.70 


-1.15 


0.85 


PW-PP 


-12.56 


0.55 


-10.25 


-6.70 


-2.58 


1.43 


-10.95 


-6.52 


-1.09 


1.02 


LDA-PAO |59| 


-12.38 


1.03 


-9.85 


-6.72 


-2.66 


1.59 


-10.63 


-6.53 


-1.14 


1.28 


EXP [0 


-13.10 


1.63 


-10.75 


-6.70 


-2.80 


2.18 


-11.24 


-6.70 


-1.30 


1.85 



TABLE III: Kohn-Sham eigenvalues calculated using various methods. The energies are calculated with respect to the top of 
the valence band and all the units are eV. PAOcxp and PAOthe are the results of the present calculation assuming the lattice 
constant to be respectively the experimental ao = 5.65A and the theoretical ao =5.635A. 



Material 


Ga (jel) 


Mn (|e|) 


As(|e|) 


As^^- del) 


GaAs 


+0.056 




-0.056 




MnAs 




-0.322 


+0.332 




Gao.938Mno.o62As 


-h0.042 


-0.089 


-0.046 


+0.005 


Gao.958Mno.o42As 


+0.046 


-0.085 


-0.049 


+0.005 


Gao.969Mno.o3i As 


+0.047 


-0.083 


-0.049 


+0.005 


Gao.979Mno.o2iAs 


+0.049 


-0.086 


-0.050 


+0.005 



TABLE IV: MiiUiken atomic charges for GaAs, MnAs and 
(Ga,Mn)As at different Mn concentrations. The last two 
columns correspond respectively to the average over the As 
atom excluding the ones coordinated with Mn and the aver- 
age over the four As atoms coordinated with Mn. The lattice 
spacing of MnAs is assumed to be ao=5.65A. 



Material 


Mn-dt (|e|) 


Mn-di (|e|) 


As-pt (|e|) 


As-Pi (le|) 


As^^-PT del) 


As"'-Pi del) 


MnAs 


4.642 


0.855 


1.365 


1.650 






Gao.938Mno.062As 


4.665 


0.788 


1.626 


1.637 


1.580 


1.647 


Gao.958Mno.042As 


4.679 


0.770 


1.628 


1.638 


1.583 


1.644 


Gao.969Mno.o3iAs 


4.675 


0.771 


1.630 


1.637 


1.584 


1.644 


Gao.979Mno.o2iAs 


4.682 


0.768 


1.632 


1.636 


1.584 


1.644 



TABLE V: MiiUiken atomic orbital populations in MnAs and (Ga,Mn)As at different Mn concentrations. The symbols | and 
I correspond to majority and minority spin respectively. The last two columns correspond to the four As atoms coordinated 
with Mn. The lattice spacing of MnAs is assumed to be ao=5.65A. 



X 


Na (eV) 


Nf3 (eV) 


1.0 


0.176 


-2.44 


0.06250 


0.272 


-5.48 


0.04166 


0.298 


-6.54 


0.03125 


0.168 


-7.34 


0.02084 


0.234 


-8.16 



TABLE VI: Exchange constants as a function of the Mn con- 
centration X for Gai-^Muj^As. 



